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Abstract 

We deal with the problem of estimating the volume of inclusions using a finite number 
of boundary measurements in electrical impedance tomography. We derive upper and lower 
bounds on the volume fractions of inclusions, or more generally two phase mixtures, using 
two boundary measurements in two dimensions. These bounds are optimal in the sense 
that they are attained by certain configurations with some boundary data. We derive the 
bounds using the translation method which uses classical variational principles with a null 
Lagrangian. We then obtain necessary conditions for the bounds to be attained and prove 
that these bounds are attained by inclusions inside which the field is uniform. When special 
boundary conditions are imposed the bounds reduce to those obtained by Milton and these 
in turn are shown here to reduce to those of Capdeboscq-Vogelius in the limit when the 
volume fraction tends to zero. The bounds of this paper, and those of Milton, work for 
inclusions of arbitrary volume fractions. We then perform some numerical experiments to 
demonstrate how good these bounds are. 
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1 Introduction 



One of the central problems of the theory and practice of electrical impedance tomography is 
the problem of estimating the volume of the inclusions in terms of boundary measurements, 
either voltage measurements when currents are applied around the boundary of the body or 
current measurements when voltages are applied. The problem can described in rigorous terms 
as follows: Let D be an inclusion inside a body fi, and suppose that the conductivities of D 
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and Q \ D are ai and (72 (cxi 7^ (72), respectively. Let a = aix{D) + (72xi^ \ ^) where x{D) is 
the characteristic function of D and the potential V be the solution to 



V • crVy = in 17, , 

V = V^ on dVL ^ ' 

for some Dirichlet data (voltage) on 517. Then the measurement of current (the Neumann 
data) is q := on d^l. (Throughout this paper ^ denotes the normal derivative.) The 
problem is to estimate the volume \D\ of the inclusion using the boundary data {V^,q) for 
finitely many voltages, say = Vi, . . . , V^. If the Neumann boundary condition = q is 
prescribed on instead of the Dirichlet condition, then the measurement is := V\dQ. 

The purpose of this paper is to consider this problem and derive optimal upper and lower 
bounds for the volume fraction of inclusions in two dimensions. In fact, we deal with a more 
general situation where 17 is a two phase mixture in which the phase 1 has conductivity ai and 
the phase 2 has conductivity (T2 (fxi > 172) so that the conductivity distribution o" of O is given 
by (7(x) = (7iXi(x) + cr2X2(x) where Xj is the characteristic function of phase j for j = 1, 2, i.e., 

f \ 1 f \ / 1 in phase 1, . . 

Xi(x) = 1-X2(x) = |q in phase 2. ^''^^ 

We derive optimal upper and lower bounds for the volume fraction /i of phase 1 (/i = 
PI using boundary measurements corresponding to either a pair of Dirichlet data 

{Vi and V2) or a pair of Neumann data {qi and 52) on 17. The bounds are optimal in the 
sense that they are attained by some inclusions or configurations. The bounds can be easily 
computed from the boundary measurements. In fact, they are given by two quantities: the 
measurement (or response) matrix A = {aij)ij=i^2 where 
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and 



^ Vil^ (1-3) 
an 



if the Dirichlet data are used. Here and throughout this paper, ^ denotes the tangential 
derivative along 517 in the positive orientation. If the Neumann data are used, then b^, is 
replaced with 

bN:=T^ [ (?i(x)( /"%2). (1.5) 



where the xq E 917 and the last integral is on the surface (917. See Theorem 12.11 and 12. 2i 

Some significant results on the problem of estimating the volume of inclusion using boundary 
measurements are as follows. Kang-Seo-Sheen [16J, Alessandrini-Rosset p], and Alessandrini- 
Rosset-Seo [2] obtained upper and lower bounds for the volume of the inclusion. However, 
their bounds involve constants which are not easy to determine, and hence it is not possible to 
compare them with the bounds of this paper. It is worth emphasizing that these results use only 
a single measurement. Another important result on volume estimation is that of Capdeboscq- 
Vogelius [BlEj. They found, using the Lipton bounds on polarization tensors [H], upper and 
lower estimates for the volume of inclusions occupying a low volume fraction, which are optimal 
bounds in the asymptotic limit as the volume fraction tends to zero. Recently it was recognised 
by Milton [25] that bounds on the response of two-phase periodic composites could be easily 
used to bound the multi-measurement response of two-phase bodies when special boundary 
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conditions are imposed (see (|2.73p and (j2.77p below) and that these could be used in an inverse 
fashion to bound the volume fraction. As shown here those bounds coincide exactly with the 
Capdeboscq-Vogelius bounds in the asymptotic limit as the volume fraction tends to zero. 

The bounds obtained in this paper allow for more general boundary conditions and we 
emphasize that they are optimal for any volume fraction. They reduce to those of Milton for 
the special boundary conditions, but have the advantage of being able to utilize the same set of 
measurements for both the upper and lower volume fraction bounds. We derive the bounds using 
the translation method which in its simplest form is based on classical variational principles 
with null Lagrangians added, i.e., non-linear functions of fields which may be integrated by 
parts and expressed in terms of boundary measurements. The translation method, developed 
by Murat and Tartar [28^ l29l |26] and independently by Lurie and Cherkaev [21^ I22j . is a 
powerful method for deriving bounds on effective tensors of composites. As shown by Murat 
and Tartar it can be extended using the method of compensated compactness to allow for 
functions more general than null Lagrangians, namely quasiconvex functions. It is reviewed in 
the books [9l [23l [3l [30]. The use of classical variational principles to determine information 
about the conductivity distribution inside a body from electrical impedance measurements was 
pioneered by Kohn and Berryman |17j . 

We continue our investigation by looking for necessary and sufficient conditions for the 
bounds to be attained. These are the exact analogs of the condition found by Grabovsky [11] 
for attainability of the translation bounds for composites. (See also section 25.6 of [23j-)- It 
turns out that the upper bound is attained if and only if the field in phase 1 is uniform and the 
lower bound is attained if and only if the field in phase 2 is uniform. It means that if phase 1 is 
an inclusion, the upper bound is attained if the field inside the inclusion is uniform. However, 
the lower bound can only be approached since no boundary data generate a nonzero uniform 
field outside the inclusion. The lower bound (for /i) can be attained for the configuration where 
phase 2 is an inclusion. 

There are plenty of inclusions inside which the field is uniform for some boundary conditions. 
We call such inclusions Efj-inclusions. They include E-inclusions which were named in [20j . An 
inclusion E is called an E-inclusion if the field inside E is uniform for any uniform loading at 
infinity. More precisely, E-inclusions are such that if V is the solution to 



then — Vy is constant in E for any direction a. If an E-inclusion E is simply connected, then E 
must be an ellipse (an ellipsoid in three dimensions). This was known as Eshelby's conjecture 
|10] and resolved by Sendeckyj in two dimensions [27] (see also |14ill9j for different proofs), and 
by Kang-Milton [15] and Liu [19] in three dimensions. There are E-inclusions with multiple 
components [H [T9l [T3] . There are also inclusions other than E-inclusions inside which the field 
is uniform. For example, if 17 contains a connected component, say E, of an E-inclusion with 
multiple components, then E is an E^-inclusion. More generally if E is an E^-inclusion and 
^' C then the field in E ^ will be uniform when appropriate boundary conditions are 
imposed at the boundary of 

We perform some numerical experiments to demonstrate how good the bounds are for 
inclusions. Special attention is paid to the variation of the bounds when certain parameters, 
such as conductivity, the volume fraction and the distance from the boundary, vary. We also 
look at the role of boundary data. 

This paper is organized as follows. In the next section we derive the lower and upper bounds 
on the volume fraction. In section 3, we obtain conditions for these bounds to be attained, and 
then in section 4, we show that if the field is uniform in phase 1 then the upper bound is 



{ 



V • {(Jix{E) + a2x(K2 \ E))VV = in 
y(x) — a • X = 0(|x|^-^) as |x| — )• oo. 



(1.6) 
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attained and if the field is uniform in phase 2 then the lower bound is attained. In Section 5, 
we obtain different sufficient conditions for the bounds to be obtained. Section 6 is devoted 
to the asymptotic analysis of the bounds when the volume fraction tends to zero. Numerical 
results are presented in section 7. In section 8 we show how to construct a wide variety of 
simply connected Ef^-inclusions, following the approach outlined in section 23.9 of 123] ■ 

We emphasize that the method of this paper (the translation method) works for three 
dimensions as well. The results in three dimensions will be presented in a forthcoming paper. 



2 Translation bounds in two dimensions 

In this section we derive upper and lower bounds on fi (the volume fraction of the phase with 
higher conductivity) using pairs of Cauchy data. Each bound requires two pairs of Cauchy data. 
The derivation in this section is based on the translation method, and parallels the treatment 
given by Murat and Tartar [28l [Ml [26] and Lurie and Cherkaev [2T][22]. 



2.1 Lower bound 

Consider two potentials satisfying 

V-aWj = inn, j = l,2. 

Let 



(2.1) 



jj(x) = -a(x)VV,(x), j = l,2. (2.2) 

We want to use information about two pairs of Cauchy data {V^,qi := — ji • n) and {V2,q2 '■= 
— j2 • n) on dQ to generate a lower bound on /i. 
Using the boundary data we can compute 



an 



(2.3) 



We assume that (ji) and (j2) are linearly independent. Then, by taking linear combinations of 
the old potentials if necessary we may assume 



(jl) 



'1 




"0" 





> (j2> = 


1 



With 



let us introduce a 4 x 4 matrix 



R± 



1 
-1 



-cR\ cj-i 



(2.4) 

(2.5) 
(2.6) 



where the constant c is chosen so that Lc{x) > for all x. Here we assume that a is an 
anisotropic conductivity (matrix). With the constants ki, k2, k^, k4^, define a 4-dimensional 
vector J by 

^iji + k2i2 
ksh + k4i2 



J 



We then consider 



(2.7) 
(2.8) 
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Define a 2 x 2 matrix A = {aij), which we call the response (or measurement) matrix, by 

(2.9) 



and 



Since 



/ ji • RUi = 0, ^ = 1,2, 



one can see that 







'k{ 


k2 




k2 


k3 


k3 






/S4 



(2.10) 
(2.11) 

(2.12) 



(2.13) 



where 

ail 112 cb 
ai2 022 —cb 
—cb ail oi2 
cb ai2 a22_ 

We emphasize that Wc can be computed from the boundary measurements. In fact, since 
V X R±ji = 0, there are potentials (pi such that 

R±3i = y^i- (2.14) 

Moreover, if t is the unit tangent vector field on dO, in the positive orientation, then 

t • = R^t • ji = -ji • n = on 50 (2.15) 

(T for the transpose) , and hence the boundary value of ipi which we denote by (p^ is given by 

<^0(x)= /%, (2.16) 

where the integration is along in the positive orientation (counterclockwise). Hence 

b = -^ I Q^^^ = ^ I a9(^V. (2.17) 



T7T7 / <ll'P2 = 77^ / (12^1- 

l"l Jan l^'l Jon 



Since 



= 7^7 / (^iji + A;2j2)cr ^(fciji + A;2j2) + (A;3ji + k4i2)cr ^{ksh + A;4j2) 
l^^l Jn 

+ 2c{kik4- k2k3)b, (2.18) 
we have the variational principle 



1 f 

mm — / 

•ji=V-j2=0 l^^lin 



V 

ji •n = gi, j^-n = q2 



kih + A;2j2 
ksh + k4j2 



(2.19) 
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One can easily see from the constraints that 

Hi) = jkl I -^1^ = Hi)- 

~ l^'l Jdn 

So if we replace the constraints by the weaker constraint that 

(ji) = (ji>, ^ = 1,2, 

then we get 



Wc > mill 

jl) j2 



ksjl + ^4j2 



kiji + k-ih 
km + ^4)2 



In order to find the minimum, we first observe that at the minimum 



kitpi + k2ilJ2 
k^Tpi + kiii2 



kiji + k2h 
ksii + ^4)2 







for any (vector- valued) functions ^i,V'2 satisfying {ipi) = (^^2) = 0, which implies 

i^c(x) 



^ijl + ^2j2 

k-sii + ^4)2 



We then have 



^iJi + ^2j2 
ksh + kij2 

Thus the minimum is given by 



fi (a constant vector). 



kih_ + k232 
ksJi + kij2 



kih_ + k2j2 
ksii + fe4j2 



which implies, thanks to (j2.4p . that 



^1 jl + ^2 j2 
fcsji + k4j2 

kiii + k2h 
kz]i + fey 2 



feljl + fe2j2 
fesjl + fe4j2 





'k{ 




'k{ 




k2 
ks 
ki 




k2 
fes 

^4 



Thus we have 

Let us now assume that a is isotropic so that 

o-^i c ■ 
cr-i -c 
- -c a-' 

c fT-i 



and 



{L 



((7-2 - c2) 



(7-1 








— c 





C7-1 


c 








c 







— c 








(7-1 
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Since 



■ 

Q'^ 



1\-1 



Q 0' 
Q 



i\-i 







for any rotation Q, we obtain from (|2.28p that 
In particular, we may choose Q so that 



Q 
Q 



an ai2 

ai2 022 



Q 



Ai 
A2 



i\-i 



(2.31) 



(2.32) 



(2.33) 



where Ai > A2 are eigenvalues of the response matrix {aij). Then by taking the inverse of both 
sides of (|2.32p we get 

r A2 -cb 
1 \i cb 

(AiA2-c262) c6 A2 
-cb Ai 



(2.34) 



So we get the inequality 
1 



(A1A2 - c262) 



A2 —cb 
—cb Ai 



V < 



[CF- 



C2) 



-c a 



(2.35) 



for any vector v. 

Now suppose that the medium is 2-phase, with ai > o"2. In this case Lc(x) > as long as 
c < a^^. We take the limit as c approaches a^^. Then 



(ar^-c2) 



0"i 



-C fJ 



1 



(2.36) 



becomes infinite unless v is proportional to 



(c7-2 - c2) 



-c a 



1 

and when v 
_ 2{a-^ - c) 

C7~2 — c2 



a^' + c 



(2.37) 



approaches cJi in phase 1 and 2/{a^ ^ + ^) in phase 2. Hence the bound in (j2.35p reduces to 

2/2 



Ai + A2 - 2b/ Gi 

2~ - ■/I'^l + 



A1A2 - 62/^2 



1/cJl + I/CJ2 
2/2'7l'72 
CTl + (T2 



_ f C^l(<7l - CT2) 2(Tl(T2 

— Jl— ^ ^ ^ H ^ ; 

(cJl + (T2j CTl + £72 

which gives the desired lower bound on the volume fraction: 



or 



/i> 
/i> 



(ai + 0-2) 
CJi(cri - (T2) 

(q-i + 0-2) 



Ai + A2-25/fTi 2(Ticr2 



A1A2- 62/(7^ CTl+(T2 

TrA-26/(Ti 2aio-2 
det A -62/0-2 CJl +0-2 



(2.38) 
(2.39) 
(2.40) 

(2.41) 
(2.42) 
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where the matrix A is defined by (j2.9p . We emphasize that the righthand side of (j2.42p can be 
computed by the boundary measurements. In fact, A is computed by using (|2.9p and b using 
(I^TT]) under the condition (f^ . 

In general, if Neumann data qi and (72 do not satisfy (|2.4p . then let 



Pn 



-1 

W\ 



an I'^l 



en 



Then ji and j2 defined by 



satisfy Since 



and 



j« = X] [^-^lim jm' i = l,2 



m=l 



1 



1 



jl • ^±j2 



detPN 



PnAPJ, 



jl • R±h 



(2.43) 

(2.44) 

(2.45) 
(2.46) 



we obtain the following theorem from ()2.42p . 
Theorem 2.1 Let P^ he given by \2.4'J^ and 

bN ■■ 

Then, 



fi> 



[ai + (72, 



X[h-R±h = TL[ qi{^){ [ 92). 
l"l JQ l"l Jan Jxo 

Tv{PnAP^) - 2(det PN)bN/(Ti 2aia2 



CJl((Tl - (J2) 



{detPNf {det A- b%/al) ai + a2 



(2.47) 



(2.48) 



2.2 Upper bound 



We now derive the upper bound on fi. 
Let us introduce a 4 x 4 matrix 



a cR± 
-cR± a 



(2.49) 



where the constant c is chosen so that L'^(x) > for all x. With the constants ki, k2, k^, k^ 
and 

e,(x) = -Vy,(x), j = l,2, (2.50) 
define a 4-dimensional vector E by 



E :-- 



kiBi + A;2e2 
k-^ei + A;4e2 



We then consider 



Ty^= {E-L',E). 
The minimization problem in this case is 



W' > min 
61,62 

(ei) = (6i) 



fciei + ^262 
A;36i + /C462 



■L'. 



kiBi + A;2e2 
A;3ei + k4e2 



(2.511 



(2.52) 



(2.53) 
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As for (j2.24p . one can show that at the minimum of the right hand side of (|2.53p 



and the minimum is given by 



fciei + /c2e2 
^361 + kie2 



H (a constant vector) 



(2.54) 



fciei + k2e2 

ksBi + /C4e2 



■L' 



kiBi + k2e2 
ksBi + k4e2 



k^ei + ^462 



UK) 



ki^i + k2e2 
/caei + A;4e2 



(2.55) 

Proceeding in the exactly same way as in the previous subsection (with c approaching to 
(T2), we can derive 'dual bounds': 



Tr^-26'a2 ^ h ^ 2/i 



where 
and 



det A — 6'2(t| CJ2 cTi + a2 
b' := (ei • R^e2) 

_ Oil a\2 
\a\2 a22 
in which 

aij ■= {ei ■ aej) = (e^ • jj), 
and hnear combination of potentials have been chosen so that 



^^^^ = W\ 



an 



on 



Apart from this constraint, ei and 62 are any fields solving 



V • aWj = inn, ej 



One can obtain from (j2.56p the upper bound on fi: 



0-2 (O-I + (T2) 



1 TvA- 2h'a2 



h' 



02 det^-6'V2 
e bounda 

^- [ V,'n.R^e2 = -j^ [ Vrt-e2 = j^[ V^- 
l^^l JdQ I'^l Jan JdQ 



(2.56) 

(2.57) 
(2.58) 
(2.59) 

(2.60) 
(2.61) 
(2.62) 



an 



(cri - 0-2) 

We emphasize that A and b' can be computed from the boundary measurements: 

1 

~ W\ 

and 

vi II • ii_L^2 — ~T7=rr / fi t • t;2 — 77=77 / ~ 

Ian 

More generally, if and V2 do not satisfy (|2.60p . then we have the following theorem in 
the same way as before. 

Theorem 2.2 Let 



(2.63) 
(2.64) 



-1 

W\ 



an I'^l 



an 



V2'n 



and 



Then 



bo 



\n\ 



an 



dt 



fi< 



'72 (O"! + 0-2) 

(ai - 0-2) 



1 Tr(P^AP^) - 2(det PD)&DfT2 



CJ2 {det Pd)^ {det A -bla^^ 



(2.65) 
(2.66) 
(2.67) 
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2.3 Special boundary data 

In the special case where the Neumann data are given by 



gi = -nji = -n 



, g2 = -nj2 = -n- 



we have 

where crj\f is the Neumann tensor which is defined via the relation 



6=1 and A = cr^^ 



(e) = cr^^(j), 



(2.68) 
(2.69) 
(2.70) 



when the Neumann data is given hy q = — n • v for some constant vector v. In fact, we have 
from (pJil) and (fZTTl) 



l"l JdQ I'^l Jan 



Ri 



1, 



and from ()2.9p 

Jan Jafi 
So, the bound (|2.42|) reduces to the bound 



(2.71) 



(2.72) 



fi> 



(di + 0-2) 
ai{ai - 02) 



Tr CTjy-'^ - 2/cri 2a\a2 



det (T^^ - 1/crf CTl + (72 



(2.73) 



of Milton [25]. 

If the Dirichlet data take the special affine form 



X, 



X, 



one can prove in the same way that 

6' = 1 and A = ctd 
where (Tij is the Dirichlet tensor which is defined via the relation 

o-D(e) = (j) 



(2.74) 



(2.75) 



(2.76) 



when the Dirichlet data V'^ is given by — v • x for some constant vector v. Thus the bound 
(I2.62P reduces to the other bound 



/i< 



02{p\ + 02) 
(cTi - 02) 



1 TV (T£) — 2(72 



(72 det ctd — a\ 



(2.77) 



of Milton [25] 
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3 Attainability conditions of the bounds 



In this section we derive conditions on the fields for the bounds in (j2.42p and ()2.62p to be 
attained. We will show in the next section that the bounds are actually attained by certain 
inclusions. 

The derivation of the lower bound on /i, and in particular (j2.24p and (|2.25p . suggests that 
if there is no column vector K = {ki, k2, k^, k^)'^, with say |Kp = kf + k2 + k"^ + kj = 1, such 
that 



L -1 



^iji + 



(3.1) 



then the lower bound will not be attained. Here, ((L^-i) ^) ^ is understood as the limit of 
{L~^)~^ as c tends to a^^. To prove this, fix cq < cr^^ and let 



F,(x) := Lc(x) 



hh + k2h 
ksh + kij2 



kih + k2h 
^sji + kih 



(3.2) 



for c such that cq < c < o^"^ . Then, we have 
(Fe • L-^F,) < (F, • L-^F,) 



CO 

kiM + k2h 
k-ih + kik 



kih + k2h 
A:3ji + kih 



ki'ii + k2h 
k-ih + k4^2 



kiji + k2j2 
ksji + kih 



(3.3) 



Letting c — t- a^^ we see that if F^-i is non-zero (in the norm) for all K with |K| = 1 then 
the right hand side of (|3.3p is non-zero in this limit, or equivalently D^-i > a((L^-i )~^)~^ for 
some a > 1. It follows that equality is not achieved in ()2.35p and hence in ()2.42p . i.e., the lower 
bound on the volume fraction is not attained. 

Conversely, suppose we have equality in (j3.ip for some K 7^ 0. Then, 



K-D -iK = K-((L -1' 



K, 



(3.4) 



and as D^-i > ((L^-i) ^) ^ it follows that D^-i — ((L^-i) ^) ^ must have zero determinant. 
A simple calculation shows that 



R± I 



where 



Hence the matrix 



'7l(o"i — (T2) 2cri(T2 



g ■= fi 



(7\ + CJ2 



0\ + C72 



0- 





Q 0' 
Q 



Ai - 1/5 
A2 - 1/5 



(3.5) 
(3.6) 



- 1/5 

-6/^1 + 1/5 





- \lg 



-b/ai + l/g Ai - l/g 











must have zero determinant, which implies 

A1A2 - (Ai + \2)/g - b'/al + 2b/{gai) = 0. 





A2 - l/g 

(3.7) 



(3.8) 
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Thus equality holds in (j2.40p and the lower bound on /i is attained. 

In summary, the attainability condition is that for some ki, k2, and /c4, 



L 



(J) 



where 



^iji + k2h 
ksh + kih 



(3.9) 
(3.10) 



Prom p.Sp a vector U in the range of {{L^-i] 



U 



ai 

-a2 
ai 



takes the form 



(3.11) 



for some oi and 02. 

We have the following theorem. 

Theorem 3.1 The attainability condition 113. 9\) for the lower bound holds if and only if 



L -iJ = U 



(3.12) 



for some U of the form 1^3. 



Proof: The 'only if part is trivial. Suppose that ()3.12p holds. We write L = L^-i for the 
ease of notation. One can see from the definition (12. 6p of L that Li (=L on phase 1) and L2 
{=L on phase 2) can be simultaneously diagonalizable. Thus in that basis (|3.12p reads as 



aS^\i(x) + A^^\2(x) j(^-)(x) = C/(^-), j = 1,2,3,4. 



(3.13) 



Here A^"*^ and Ag"''' are eigenvalues of Li and L2, respectively, and J^^\x.) and U^^^ are j-th 
components of J and U in new basis. 

Since Li has rank 2, two of eigenvalues A^''^ are zero, say A^^'' and and hence xiJ^^H'^) 
and X2J^^\^) may depend on x. However, J^^\x) for j = 1,2 is piecewise constant, and by 

(EH, 



J(^')(x) 



J ?7(J) /X[^^ in phase 1 
\ /A^-''^ in phase 2. 



Thus we have 



(j(^-))= /i/aS^V/2/A^^'^ t/(^-) = (L-i),,C/W, i = l,2. 



,(i) 



Here {L is the (j,j)-entry of the diagonal matrix (L ^). So, 

= ((L-i)-i),,(j(^-)), J = 1,2. 



(3.14) 
(3.15) 

(3.16) 



If j = 1,2, then ((L ^) ^)jj = and U, which belongs to the range of {L ^) ^ satisfies 
U^^^ = C/(^) = 0, and hence (|3.16p holds for all j. Therefore 

U = (L-i)-i(J), (3.17) 

and hence (13.91) holds. I 
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Similarly one can show that the attainability condition for the upper bound is that for some 
ki, k2, and k^ 

L' E=((L' )-i)-i(E) (3.18) 



where 



and it is equivalent to 



E 



kiGi + k2G2 

k-iBi + /c4e2 



(3.19) 



(3.20) 



for some U of the form (j3.1ip . We emphasize that the attainability conditions p.l2p and (j3.20p 
are precisely analogous to those found by Grabovsky [11] for composites. 



4 Attainability and uniformity 

We now investigate the attainability condition more closely. (|3.2Up says that the field E is 
uniform in phase 1. This condition alone guarantees that the upper bound is attained. In fact, 
we show in this section that even more is true: if the field is uniform in phase 1 for a single 
boundary data = Vi then there is a V2 such that the upper bound is attained. 

Theorem 4.1 Suppose that phase 1 and 2 have finitely many connected (possibly multiply 
connected) components and the interfaces are Lipschitz continuous. Let V be the solution to 



V • aW = inn 

V = V° on dn. 



..^n (4-1) 



in phase 1 for some 



// — VV^ is constant (the field is uniform) in phase 1 for some boundary data = 7^ 0, 
then there is a V2 such the upper bound is attained. 

Proof: Phase 1 can be broken into connected components a = 1, 2, . . . , m, and phase 
2 can be broken into connected components ^'2'^^ /3 = 1, 2, . . . , n. If ^'2^^ has a boundary in 
common with ^^"^ we denote the common boundary by F"^. 

Let Vs(x) denote the potential V{x) inside ^'f \ If -VV = 

constants ei and 62, then 

y(x) = -eix - e2y + c„ (4.2) 

for some constant Cq, inside "if^^ (where Cq = if ^'["^ touches ^'[''^^ at a common point), and 
the continuity of the potential on F"^ implies 

Vp{x) = -eix - e2y + ca on F"^. (4.3) 

Since V • j(x) = in 0, there is a continuous potential W^(x) such that 

j(x) = -CT2i?±Vl^(x) in n. (4.4) 



In phase 1, inside we have 



VVF(x) = — 

0-2 



62 

-ei 



(4.5) 



and hence 

W{^) = —{e2X-eiy) + dc, (4.6) 

0"2 
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for some constant da (where, by continuity of the potential W, da = d^ if ^^"^^ touches "^^^^ at 
a common point) 

Let VF^(x) denote the potential l^(x) inside ^'2'^'*. Since W{x) is continuous, 

'^^-{e2X - eiy) + da on T"'^. 



W-^(x) 



0-2 



Note that inside 5' 



2 ' 



0-2 



■j(x) = R^VW{^), 



(4.7) 



(4.8) 



i-e., Vp^x = WjS^y and Vp^y = —Wp^x-, which are the Cauchy-Riemann equations. Thus Vp + iW^ 
is an analytic function oi z = x + iy. 
Now consider 



-Wf3{^) + {ai/a2 + l){e2X-eiy) (4.9) 

Ty^(x) := F^(x) + (ai/a2 + l)(eix + ezy). (4.10) 

Clearly 

+ iW'p = i{Vp + iW^^) + (c7i/a2 + l)(e2 + iei){x + iy) (4.11) 
is an analytic function of z. On F"'', we have 

= -{ai/a2){e2X - eiy) + da + icri/a2 + l)(e2X - eiy) = e2X - eiy + da, (4.12) 

W'^ = -eiX - 62?/ + Ca + ((7i/c72 + l)(eiX + 62?/) = ((7l/(72) (ciX + 62?/) + Cq. (4.13) 

So the conductivity equation V • aW = is satisfied with potentials V and W defined by 
V' = V^, W = VF^ in ^' 



/9 



■2 ' 



and 



in ^^^^ . Note that 



F'(x) = 62^ - eiy + da, PF'(x) = {ai/a2){eix + e2y) + Ca 
-VV' 



-62 

ei 



(4.14) 
(4.15) 



We then have, in ^2 



(/3) 



(72/ CJ2-R_L 




'-VV' 




(T2i?± (72I 




-W 





02! (T2-R_L 
— (T2-R± (72/ 



VW;3 - ((7i/cr2 + 1) 



62 

-ei 



-0-2 (VVg - i2±VW^) + (c7i + 0-2 
a2(i?±Vy;3 + VW/3) + (cJi + (T2) 



ei 






ei 


-62 




= (0"1 + CT2) 


62 

-62 


^1 . 









(4.16) 



and in phase 1 



ail a2R± 
-a2R± (7il 



'-VV' 




_-VV'_ 





ail a2R± 
-a2R± ail 



61 
62 
-62 

61 



((71 + (72 ) 



61 
62 
-62 
61 



(4.17) 
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Thus ~ ^ where U is of the form (j3.1ip . Hence the upper bound is attained when we 

take boundary data V^^ = V'^ and = V"^ I 

Observe that the Dirichlet condition in (|4.1|) may be replaced with the Neumann condition. 
One can prove in the exactly same way that the lower bound is attained if the field is uniform 
in phase 2. 



5 Attainability and analyticity 

We have seen that uniformity and independence of the fields ei = — VVi and e2 = — VV2 in 
phase 1 is necessary and sufficient to ensure that the upper bound is attained. Now we will 
see there is a condition on the potentials Vi and V2 in phase 2 which is also necessary and 
sufficient to ensure that the upper bound is attained. We assume that phase 2 is connected and 
completely surrounds each inclusion of phase 1. 

First suppose that the upper bound is attained. Then, given a constant k, there exist 
potentials V and V', which are linear combinations of the potentials Vi and V2, such that in 



phase 1 —W 



and -VV 



Thus the analysis of the previous section holds with 





k 

ei = k and 62 = 0. In particular, we may choose k = l/(cri/cr2 + 1) and, since in phase two 

V + iW is an analytic function of z, it follows from (14. 9p that V — iV' + x is an analytic function 
of z = X + iy in phase 2. 

Conversely suppose there exist potentials V and V' , which are linear combinations of the 
potentials Vi and V2, such that V — iV + x is an analytic function of z in phase 2. Then the 
harmonic conjugate to V in phase 2 \s —V' — y and the harmonic conjugate to V' in phase 2 is 

V + X. Since by ()4.8p these harmonic conjugates can be identified with the potentials W and 
W , we have in phase 2 

W = -V'-y, W' = V + x, (5.1) 

and in particular these identities hold on the boundary of an inclusion of phase 1. By (j4.4p 
inside that inclusion V + i{(T2 / cri)W and V' + i{a2/cri)W' are analytic functions of z. Therefore 

V + iia2/ai)W' - iiai/a2)iV + i(a2M)PF) - i(x + iy) 

is also an analytic function of z inside the inclusion and from ()5.ip takes the value 

i{a2/(Ji - l)[{ai/a2 + 1)V + x] 

at the boundary of the inclusion. Since the only function which has zero real part at the 
boundary of a closed curve, and which is analytic in the interior, is an imaginary constant, 
we deduce that (o"i/(T2 + 1)^ + x is constant around the boundary of the inclusion and hence 

Ik] 

constant inside, i.e., in the inclusion —W = ^ with k = 1/((Ti/o"2 + 1)- 

The harmonic conjugate to V inside the inclusion is then —ky which can be identified with 
{cr2/cri)W (to within an additive constant). Then from the first condition in (j5.ip it follows 
that (to within an additive constant) V takes the value —ky around the boundary of the 

"0" 



inclusion and hence in its interior too, i.e., in the inclusion —W 



k 



Hence the uniform 



field attainability condition is met, and the upper bound is attained. 
We summarize our findings as a theorem. 

Theorem 5.1 Provided the body Q consists of inclusions of phase 1 completely surrounded by 
phase 2 then the upper bound is attained if and only if there exist potentials V and V , which 
are linear combinations of the potentials Vi and V2, such that V — iV + x is an analytic function 
of z = X + iy in phase 2. 



15 



Similarly we have the following theorem for the lower bound. 

Theorem 5.2 Provided the body 0, consists of inclusions of phase 2 completely surrounded by 
phase 1 then the lower bound is attained if and only if there exist potentials V and V' , which are 
linear combinations of the potentials Vi and V2, such that V — iV + x is an analytic function 
of z = X + iy in phase 1. 



6 Asymptotic bounds for small volume fraction 

Suppose that the phase 1 occupies a region uj C 0, satisfying 

dist{uj,dn) >c (6.1) 

for some c > 0. The purpose of this section is to compare the bounds (|2.42p and (j2.62p with 
the bounds obtained in [7] when the volume |ti;| of uj tends to 0. 

Let g be a function on L'^(dCl) satisfying Jq^ q = 0. Let V be the solution to 



V • aVV = in 0, 

(6.2) 

V = UJ, 

'an 



dV f 
a—=q ondn, (/ V = 0), 
9n Jen 



and let U be the solution to (j6.2p with a replaced with cj2. It is proved in [6] that given a 
sequence Un satisfying (j6.ip and such that |a;„| — )■ there is a subsequence still denoted a;„, a 
probability measure d^ supported in the set {x \ dist(x, dO.) > c}, and a (pointwise) polarization 
tensor field M(x) such that if is the solution to (j6.2p when u = cvn, then 

K(x)-C/(x) = -|u;„| / Vf/(z)-M(z)V^iV(x,z)(i//(z)+o(|u;n|), x G 90, (6.3) 
Jn 

where N{x,z) is the Neumann function for ri, i.e., U is given by 

U{z)= [ N{x,z)q{x)ds{x). (6.4) 
JdQ 

Note that we have absorbed a factor of cJi — (T2 into the definition of M given by Capdeboscq 
and Vogelius to be consistent with the conventional definition of polarization tensors. 
Let be the solution to 

' V • aVV = inn, 

v = v^ ondn ^^-^^ 

with UJ = iOn and U' be the solution to (16. Sp with a replaced with <T2. Then we have 
8V' dW f 8 

V^(x) - CT2 — (x) = \ojn\ / VU'iz) ■ M(z)V,^G(x, z)d/i(z) + o(|a;„|), x E OQ, (6.6) 
on dn Jq dn^ 

where G(x, z) is the Green function for 0, i.e., U' is given by 

U'{z)= [ J-G{x,z)V'{x)ds{x). (6.7) 
Jan 

To see ()6.6p let us define the Neumann-to-Dirichlet (NtD) map by 

Aa[q] := Vlan (6.8) 
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where V the solution to (j6.2p . Let Ag-j be the NtD map when a is replaced with a2. Observe 
that because of (16.41). we have 



V[/(z) • M{z)V,N{x,z)dfi{z) 
So (|6.3p can be rewritten as 



V,iV(y,z)-M(z)V,iV(x,z)d;u(z) 



■.= K[q]i^). 



Aa[q] = A^Jg] - \uJn\K[q] + o(|tJ„|). 
Then the Dirichlet-to-Neumann map A^^ is given by 

= {I- \un\A-,'K)-^A-,' + = A^^^ + |a;„| A^^^i^A^^^ + o(|^„|). 

Observe that 



A;^i[iV(.,z)](x) = -^G(x,z), xgJ], zea 
dnx 



9(y)f^s(y) 

(6.9) 

(6.10) 

(6.11) 
(6.12) 



In fact, 



/ A;^i[iV(.,z)](x)FO(x)= / iV(x,z)a2 — (x) = C/'(z) for ah z G f), 
Jan JdQ o'n 



(6.13) 



and hence (j6.12p follows. We now obtain (j6.6p from (j6.1ip . 

Let C/i(x) = — (T^ X and C/2(x) = — o"^ and let Vj be the solution to (j6.2p with g = qj for 



j = 1,2, where qj = (^2-^^ = ~f^j- Then, we have 

= aij = 7^ I Viqj = 7^ I {Vi- Ui)qj + a^^6, 
I'^l Jan I'^l Jan 



■n 



where is Kronecker's delta. Since 



N{x,z)qj{x)ds{x) = Uj{z), 



we have from (16.31) that 



1 

W\ 



an 



{Vi - Ui)qj = -\ujn 



Mij{x.)dii{x) + o{\uJn\)- 



Thus we have 
where 

We then have 



M := / Mij{x)dfi{:K.). 
Jn 

Tra-^'-2a:[\ -ff Tr M + ^ _ A + 

det tr^i - _ A Tr M + ^ - ^ + o(/i) 



(6.14) 
(6.15) 

(6.16) 

(6.17) 
(6.18) 



2(<T1-CT2) 
CTIO'2 


fl ^ /^''i ,TrM 

i(T2(CTl-CT2) 


+ o{fi) 




1 A'^' ij, Tr m] 


+ o{fi) 



2aia2 



1 + P-^ -TrM 

2(72 (CTl + 0-2) 



+ o(/i), 



(6.19) 
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and hence 



(o-i + (J2) 



det 



1/af 0-1 + (72 



-p^TrM + oih). 



The lower bound ()2.73p now reads 



or equivalently 



(6.20) 

(6.21) 
(6.22) 



(at - af, 

up to o(/i) terms by ()6.17p . 

We now consider the upper bound. Let C/i(x) = —x and ?72(x) = — y, and let Vi be the 



solution to (j6.5p with = C/j on 50 for i = 1, 2. Then, defining qj = 0'2^- on dVL, we have 



Wuh = aij = 1^ / Vl^Qj = 1^ / - ^2^) + f^2<Jij 

I I "/ I I "/ 



One can use (16.61) and the fact that 



to derive that 
Thus we obtain 



/ /-G(x,z)yP(x)dx = C/,(z) 
Jan 

(TD = flM + (J2I + o(/i) =CTN + o{h] 

1 Trcr£,-2o-2 /idetM 



CJ2 detcTD - (t| TrM 



+ o(/i). 



Since detM = (Tr M~^)~i Tr M, ^0711^ reads 



/i< ^t^+"i (TrM-^)-\ 



or equivalently 



„ (CTI + (T2) r , -1 \ 

/i < 7 Y(Tr(-^ + f^2 '^i'J 

(cJi - 0-2) 



(6.23) 

(6.24) 

(6.25) 
(6.26) 

(6.27) 
(6.28) 



up to o(/i) terms. 

By (f6T7l) and (fOSjl . we have 



modulo o(/i). Hence by putting (I6.22|) and (I6.28P together, we have 

2T Tr(/ - (T2(T^ )<fi< 7 r Tr(/ - (J20- 

(af - a^) CTi - CT2 



N 



(6.29) 



(6.30) 



modulo o(/i), where <t^^ is determined from the boundary measurements with special Neumann 
conditions, via ()6.14p . We emphasize that these asymptotic bounds for small volume fraction 
were found in [B] 1^ . From (j6.29p we also have the bounds 



2 2\ -'-^^'^2 ^D-l-)<h<-, 7(0^2 ^D-^) ) 



(a2 - al) 



Gi - (J2j 



(6.31) 
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modulo o(/i), where <T£) is obtained from the boundary measurements with special Dirichlet 
conditions. 

It is interesting to observe that the translation bounds also yield the Lipton bounds for the 
polarization tensor: We obtain from (j6.2ip and (j6.27p that 

TrM<^^^^ and T,{M-^) < ^^1±^^. (6.32) 

(Tl 0-2(cri - (T2) 

We refer to [U [5l [23] for properties of polarization tensors. If the phase 1 is an inclusion (or a 
cluster of inclusions) of the form 

D = eB + z (6.33) 

where e is a small parameter representing the diameter of D, B is a reference domain containing 
0, and z indicates the location of D inside Q,, then Mij = \B\~^M[B) (a constant matrix) and 
d/j, = lim„_>.oo \ujn\~^x{^n)d:>i. = 5{x — z)dx. Here M{B) is the polarization tensor associated 
with B. Therefore we have M = \B\^^ AI{B), and hence 

Tr(M(5)) < (6.34) 

0"! 

and the lower bound is given by 



Tr(M(S)-i) < (6.35) 

The bounds in (I6.34|) and (I6.35P were obtained by Lipton [18] and later by Capdeboscq-Vogelius 
[6l [7] in a more general setting. They can also easily be derived from the bounds of Lurie and 
Cherkaev [21] and Tartar and Murat 1 26\ |29] using the observation made by Milton [24] that 
the low volume fraction limit of bounds on effective tensors of periodic arrays of well-separated 
inclusions yield bounds on polarization tensors. We also mention that if the lower bound in 
()6.35p is attained for B and B is simply connected, then B is an ellipse. This was known as the 
Polya-Szego conjecture and resolved by Kang-Milton fl4l[l5] (see also a review paper [12]). 



7 Numerical results 
7.1 Forward solutions 

We implement an integral equation solver in FORTRAN in order to generate forward solutions 
of the Neumann and Dirichlet problems of the equation V-crVy = in when D is an inclusion 
and a = aix{D) + (^2X{^\I^)- We set (T2 = 1 throughout this section. We compute the forward 
solutions V with N = 64,80,96,120,160,192,240,320 and 480 equi-spaced points on dD and 
points on d^l. And then they are computed with the solutions on the finer discretization of 
A^ = 960. Figure [1] shows the convergence of a forward solver for the Neumann problem as a 
function of discretization points, N, and Figure [2] for the Dirichlet problem, with ai = 10 



7.2 Numerical Experiments 

We perform numerical simulations to judge the performance of the bounds when relevant pa- 
rameters are varying. Parameters under consideration are the conductivity contrast 0"i/cj2, the 
volume fraction /i, and the distance between the inclusion and dO,. We also investigate the 
role of boundary data in deriving bounds. 
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10' 10^ 10^ 



Number points per boundary 

Figure 1: Convergence error of the forward solver with 64-480 discretization points. The soUd 
hne represents the convergence error of V on for the Neumann problem. 



io' 




10' 10^ 10^ 

Number points per boundary 



Figure 2: Convergence error of the forward solver with 64-480 discretization points. The solid 
line represents the convergence error of ^ for the Dirichlet problem. 



We use boundary data of special forms; qi 



n and q2 



■ n as Neumann data 



for the lower bound, and Vi 



X and V2 



■ X as Dirichlet data for the upper 



bound, in all examples except examples 17.41 and 17.51 Thus except in these examples, the bounds 
correspond to those derived by Milton [25\ . 
Let 



+ £72 ) 



TvA-2b/ai 2aia2 



denote the lower bound on /i and let 



U{ai) 

denote the upper bound on fi. 



det A- b'^/af ai + a2 

TtA- 26' £72 



IfJl - £72 ) 



1 



£72 del A -b'^a?: 



(7.1) 



(7.2) 



Example 7.1 (variation of £7i). We compute the bounds changing £7i, keeping £72 = 1, when 
the inclusion is a disk or an ellipse inside a disk or a rectangle (with corners rounded). Figures 
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[3l HI and [5] show the numerical results. Figure [6] is when the inclusion is simply connected and 
of general shape. Figure [7] is when the inclusion is not simply connected. The results show that 
the lower bound deteriorates seriously as the conductivity ratio ai increases while the upper 
bounds are relatively good even with large ai. 



1.2 




5 10 15 20 

c 



1.4 




5 10 15 20 

o 



first diagram second diagram 



CTl 


L(<yi)/fi 


U(,Tl)/h 




CTl 


L{<yi)/fi 


U{<Tl)/h 


1.1 


0.9979 


1.0000 




1.1 


0.9904 


1.0077 


1.2 


0.9925 


1.0000 




1.2 


0.9783 


1.0149 


1.5 


0.9635 


1.0000 




1.5 


0.9340 


1.0337 


2 


0.8979 


1.0000 




2 


0.8532 


1.0583 


3 


0.7673 


1.0000 




3 


0.7115 


1.0917 


5 


0.5787 


1.0000 




5 


0.5237 


1.1287 


10 


0.3518 


1.0000 




10 


0.3113 


1.1659 


20 


0.1958 


1.0000 




20 


0.1710 


1.1889 



Figure 3: The bounds with increasing cxi when the inclusion is a disk and f2 is a circle, and 
fi = 0.09. We take Neumann and Dirichlet data of the special forms. The second and third 
columns are graphs of the same data; the third column is with a log-scale for the ai-axis. The 
values for the bounds are given in the table. 

Example 7.2 (variation of /i). We compute the bounds for various volume fractions. Figure 
[8] shows the numerical results. It clearly shows that the lower bound works better for higher 
volume fractions. 

Example 7.3 (variation of distance from d^l). We compute the lower and upper bounds 
changing the distance between the inclusion and dfl. Figure [9] shows the numerical results when 
o"! = 2. It shows that the further the inclusion is from dQ, the better bounds are. 

Example 7.4 (boundary data). In the example we compute the bounds using other bound- 
ary data. We use as Neumann data for the lower bound qi = —ni — nin2 and q2 = —n2 — nin2, 
and as Dirichlet data for the upper bound Vi = —x — xy and V2 = —y — xy. Figure [TOl shows 
that the special boundary data work much better. 
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1.5 




5 10 15 20 

a 



first diagram. second diagram 











CTl 


i(<Jl)//l 




1.1 


0.9982 


1.0000 




1.1 


0.9921 


1.0062 


1.2 


0.9934 


1.0000 




1.2 


0.9819 


1.0119 


1.5 


0.9677 


1.0000 




1.5 


0.9435 


1.0268 


2 


0.9091 


1.0001 




2 


0.8712 


1.0459 


3 


0.7895 


1.0001 




3 


0.7396 


1.0714 


5 


0.6099 


1.0001 




5 


0.5569 


1.0988 


10 


0.3818 


1.0002 




10 


0.3395 


1.1257 


20 


0.2170 


1.0002 




20 


0.1896 


1.1420 



Figure 4: The bounds with increasing cji when the inclusion is an elhpse and is a circle and 
/i = 0.08. We take Neumann and Dirichlet data of the special forms. The second and third 
columns are graphs of the same data; the third column is with a log-scale for the cJi-axis. The 
values for the bounds are given in the table. 
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first diagram second diagram 



CTl 


L{Ti)/fi 






CTl 




C/{<Ti)//i 


1.1 


0.9976 


1.0003 




1.1 


0.9915 


1.0065 


1.2 


0.9917 


1.0006 




1.2 


0.9803 


1.0125 


1.5 


0.9614 


1.0013 




1.5 


0.9376 


1.0281 


2 


0.8939 


1.0022 




2 


0.8576 


1.0480 


3 


0.7608 


1.0033 




3 


0.7155 


1.0744 


5 


0.5708 


1.0044 




5 


0.5262 


1.1027 


10 


0.3449 


1.0054 




10 


0.3122 


1.1302 


20 


0.1912 


1.0060 




20 


0.1713 


1.1468 



Figure 5: The bounds with increasing ai when the inclusion is a disk and is a square, and 
/i = 0.0699. We take Neumann and Dirichlet data of the special forms. The second and third 
columns are graphs of the same data; the third column is with a log-scale for the cJi-axis. The 
values for the bounds are given in the table. 



1.2 




5 10 15 20 



Figure 6: The bounds with increasing cji when the inclusion is not a disk or an ellipse and il. 
is a square, and /i = 0.0673. We take Neumann and Dirichlet data of the special forms. The 
second and third columns are graphs of the same data; the third column is with a log-scale for 
the (Ti-axis. 
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Figure 7: The bounds with increasing cji when the inclusion is an annulus. We take Neumann 
and Dirichlet data of the special forms. The second and third columns are graphs of the same 
data; the third column is with a log-scale for the ui-axis. 
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Figure 8: ai = 5. The bounds changing the volume fraction. We take Neumann and Dirichlet 
data with the special forms. 




Figure 9: a"i = 2 and fi = 0.0262. The bounds changing the distance between the inclusion 
and 50. We take Neumann and Dirichlet data of the special forms. 




10 20 30 



Figure 10: The bounds changing cji in case that we take Neumann data gi = —ui — v\V2 and 
52 = —1^2 — Also we take Dirichlet data V\ = —x — xy and V2 = —y — xy. 
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Example 7.5 When we use special Neumann data qi = —ni and q2 = —n2, then a pair of 
Dirichlet data are measured on d^l. We may use this data to compute the upper bound using 
the formula (|2.67p . Likewise, we may use the measured Neumann data corresponding to the 
Dirichlet data Vi = —x and V2 = —y to compute the lower bound using the formula ()2.48p . 
Figure [TT] shows numerical results when the volume fraction varies. In this example it clearly 
shows that bounds using the measured data are better than those using the given data. 




Figure 11: ai = 5. LN{fi) is the lower bound using the Neumann data corresponding to the 
special Dirichlet data and Unifi) is the upper bound using the Dirichlet data corresponding to 
the special Neumann data. 



8 Construction of Ej^-inclusions 

Following the method outlined in section 23.9 of |23) we look for a simply connected inclusion 
inside which the field is uniform for some boundary condition assigned on the outer boundary. 
More precisely, we look for an inclusion E contained in a domain (bounded or unbounded) 
such that — Vy is uniform inside E, where V is the solution to 

V-aVV = inn, 

V = V^ on dn ^ ' 

for some boundary data with a = aix{E)-\-a2x{^\E) {ai 7^ (T2). We may suppose, without 
loss of generality, that e = —W = (—1,0)^. We also suppose that the coordinates have been 
positioned and scaled so that ymax = 1 and y^[-^ = —1, where ymax = max{y | (x,y) G 
E for some x} and y^[-^ = min{y | (x, y) £ E for some x}. Let W he a harmonic conjugate of 
y in $7 \ so that V + iW is an analytic function of z = x + iy in Q\ E. Then we have 

V = x, W = — y, on dE. (8.2) 

0"2 
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Define new potentials u and v by 



u + iv = 



i{V + iW -z) 

(1 - (Ji/(J2) 



(8.3) 



Then, u + iv is still an analytic function of z = x + iy in Q \ E, and on dE 



u = 



-W + y 
1 - ai/cJ2 



= y, 



V 



1 - CTl/cJ2 



V -X 



= 0. 



(8.4) 



Now assume n + it; is a univalent function of z = x-\-iy inside Vt\E, and consider x + iy as 
an analytic function of u + iv (hodograph transformation) . Because of ()8.4p , the image of dE 
by n + if is the slit S = [yj^i^' J/max] = [—1, 1] on the n-axis, and y = u on S. 

The problem is now to construct a function z = x + iy = f{u + iv) such that 

(i) / is analytic and univalent \nU \S for some neighborhood U of S, 

(ii) Im/ = on S, 

(iii) Re /|+ — Re /|_ > on 5 except at ±1 where it is 0. 

Here |+ and |_ indicate the limit from above and below 5, respectively. 

One can see that the conditions (i) , (ii) , and (iii) guarantee that / maps U \ S onto Q\E 
for a simply connected domain E and Q a domain containing E. In fact, (ii) and (iii) imply 
that / maps S onto dE and the orientation is preserved. Since / is conformal, it maps U \ S 
to outside E. 

We have the following lemma for univalence. 

Lemma 8.1 Let j be a simple closed curve which consists of two curves 7^ and 7^. Let U be 
an open neighborhood of S and let Bi(5) and B^i{5) be open balls of radius 6 centered at w = 1 
and w = —1 respectively. Let f be an analytic function in U\S which maps U \ S outside 7 of 
the form 



where Img = on S. Suppose that the mapping u 1— )• limj,_^o+ f{u + iv) is one-to-one from S 
onto 7"*", and u 1— )• limj,^o+ /(u — iv) is one-to-one from S onto 7^. If there is 6 > such that 
f is univalent in Bi{6) \ S and in B^i(5) \ S, then there is an open neighborhood Uq of S such 
that f is univalent in Uo \ S. 

Proof: Let ip{z) = {z -\- ^)/2 for \z\ > I. ip maps \z\ > 1 onto C \ S and \z\ = 1 onto S. 
Let G{z) = g[(p{z)). Since ImG(2:) = on \z\ = 1, G can be extended so that it is analytic in 
1 — e<|^;|<l + efor some e > 0. Let F{z) = f{ip{z)). Then F is analytic \n 1 — e <\z\ < 1 + e 
and univalent in the neighborhoods of z = 1 and z = —1. Moreover, F is one-to-one from 
\z\ = 1 onto 7. We claim that F is univalent in 1 — < l-^l < 1 + foi' some £q > 0. In fact, 
if not, then for each n there are zi^n and Z2,n such that 1 — ^ < \zj^n\ < 1 + ^i,n / -^2,™, and 
F{zi^n) = F{z2,n)- For j = 1, 2, the sequence Zj^n has a subsequence which converges to a point 
on \z\ = 1, say zj. Since F is one-to-one on \z\ = 1, zi = Z2. But this implies that F'{zi) = 0, 
where 



and since g'{Lp{zi)) is real we conclude that zi = 1 or zi = —1 which is contradiction since F 
is univalent in the neighborhoods of these points. Thus F is univalent in 1 — eo < < 1 + eo 
for some eo > 0. This completes the proof. I 



f[w) = iw + g{w) 



(8.5) 



F'{z) = f'{v{z))^'{z) = [i + g'{v{z))]{l - z-')/2 



(8.6) 



27 



We now construct / satisfying (i), (ii), and (iii) using conformal mappings. Let w = u + iv 
and define 

9{w) = f{w) - iw (8.7) 

so that Im^f = on 5. Let 

which maps S onto the positive real axis. Let C, = with the branch cut along the positive 
real axis and define 

m = 9 (^) . (8.9) 

Then ImF = on the whole real axis. Thus, by defining F{C,*) = F{Q*, where * denotes the 
complex conjugate, F can be extended as an analytic function in a tubular neighborhood of 
the real axis. Moreover, since g is analytic in a neighborhood of —1 except the part of the slit 
and the bilinear transform maps a neighborhood of —1 onto outside a compact set, F must 
be analytic in C \ (K U K*) where if is a compact set in the upper half plane and K* is its 
symmetric part with respect to the real axis, i.e., K* = {z* \ z G K}. F satisfies 

(i) ' F is analytic in C \ (if U K*) for a compact set K in the upper half plane. 

(ii) ' ImF = on the real axis, 

(iii) ' F{C) - F{-C) > for real positive C- 
The function / is now given by 



f{w) = iw + F\^^^-^.] (8.10) 

Note that y = u on the slit and hence dE is given by 




In addition to (i)', (ii)', and (iii)', F needs to be univalent inside a sufhently small ball around 

the origin, and outside a sufficiently large ball. The first condition is satisfied if F'{0) ^ 0. Since 
F maps oo to a point in C, F being analytic and univalent outside a sufficiently large ball has 
the series expansion 



m = E 7J (8-12) 

j=0 ^ 



as — > oo, where /3i / (and Pi is real and positive from conditions (ii)', and (iii)') 
We make a record of these conditions: 

(iv)' The derivative F'{0) is non-zero, and F{() has the asymptotic expansion 



/3: 

= Po -I- 

where /3i is real and positive. 



i^(0=/5o + ^ + 0(|Cr') aslCl^oo, (8.13) 
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Good candidates for functions satisfying (i)', (ii)', and (iv)' are rational functions of the 
form 

" 6r, b* ' 



a=l 



+ 



C-aa C - a* 



+ c 



^.14) 



where the a^'s are complex numbers with positive imaginary parts, the b^s are complex num- 
bers, c is a real number, and 



E Re(6„) > 0, Yl Mba/al) + 0. 



^.15) 



OL = \ 



To ensure that (iii)' is satisfied we require that the function 



F{C) - F{-C) = 



+ 



b* 



.16) 



has no real roots aside from ^ = 0. (The sign of the inequality in (iii)' is guaranteed by the 
positivity of /3i.) 

Let us now characterize those rational functions F which yield ellipses as E^-inclusions. 
Because y = u on the slit [—1,1], the ellipse takes the shape like the first figure in Figure O 
(after translation). Let the ellipse be given by + ay'^ + = c with 4a > /3^. Solving for x 
we get 

_ -Py ± _ 4(a2/2 - c) 



Since the discriminant vanishes at y = ±1, we have c = 4a — /?^, and hence 

+ 



X 



2 ^ 2 



(4a - /32) 



1 + y' 



^.17) 



(8.18) 



Letting C = \/ 1+| > we have 



X 



— eT~i — + 2 =^^^) 



for real C. It means that ellipses are obtained by F's of the form 

b b* 



F{C) 



with a = i and b with positive real part. 



C-a C - a* 



+ c 



(8.19) 



?.20) 



Example. In this example, we construct some En-inclusions other than ellipses. We use F in 
the form (]8.20p with c = (it amounts to translating the figure). Then in (^-coordinates / is 
given by 

^<« = ?TT + (^ + C^- 

where both (j8.15p and the absence of real non-zero roots of (j8.16p will be ensured if we choose 
b and —bja^ with positive real parts. 

We will plot the image of a vicinity of the real axis in the upper half plane under the map /. 
To avoid computational difficulty in dealing with an infinite space, we use a bilinear transform 



c 



1 — iw 



w — I 



1.22) 
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which maps the unit disk onto the upper half plane. Then we plot 



2i 



+ 



CH2+1 CH 



+ 



b* 



C{w) 



i.23) 



for w = re*^ with 1 — e < r < 1. From the expansions for -F(C) hi powers of and 1/C we see 
that near the bottom and top of the inclusion the boundary is given by 



Re(6)V2(l + y)+0(l + 2/), 



-2Re{b/a) - Re{b/a^)^y2{l - y) + 0(1 - y). (8.24) 



Thus the bottom and top are positioned at x = and x = — 2Re(6/a) and the curvature of the 
boundary there is determined by Re(6) and — Re(6/a^) respectively. 

Figure [T3l shows various shapes of d^l, which are the image of |2;| = r < 1 under /, and the 
boundary of En-inclusion, which is the image of |2;| = 1. Figure [T ^ [T ^ [TCI [TUl [T71 and [TH] show 
various shapes of EQ-inclusions when we vary the complex parameters a, b, and 

We emphasize that with these values of a and b, the univalence of / is guaranteed by Lemma 

EE 



a=0.0+i 



a=0.4+i 





a=0.8+i 




-1 1 



Figure 12: Various shapes of E^-inclusions when varying Re a with Im a = 1 and 6=1. 



0.5< |zl=r< 1 f{z) on |zhr 




Figure 13: With a = 0.8 + i and 6=1, the inner most curve (the image of \z\ = 1) is the 
boundary of the E^-inclusion (the rightmost inclusion in Fig. [T2]) . The others are images of 
\z\ = 0.9, 0,8, 0,7, 0,6, 0,5. These, or any simple closed curve enclosed by them, can be 
regarded as boundaries of 0. 
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Figure 14: Various shapes of En-inclusions when varying Im a with Re a = and b = 1 + i. 



b=0.4+i b=1.0+i b=1.6+i 




-3 -2-1 1 -3 -2-1 1 -3 -2-1 1 

Figure 15: Various shapes of E^-inclusions when varying Re b with Im 6 = 1 and a = 1.3i. 

b=1 .0+0.4 i b=1.0+1.0i b=1.0+1.6i 




-3 -2-1 1 -3 -2-1 1 -3 -2-1 1 

Figure 16: Various shapes of E^-inclusions when varying Im b with Re 6 = 1 and a = 1.3i. 



Imb=1.0 lmb=7.0 lmb=10.0 lmb=15.0 




-40 -30 -20 -10 -40 -30 -20 -10 -40 -30 -20 -10 -40 -30 -20 -10 

Figure 17: Various shapes of En-inclusions when varying Im b with Re b = 0.1 and 6/a^ = 



-10 + 2i. 
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Imb/a^=1.0 lmb/a^=5.0 lmb/a^=10.0 lmb/a^=15.0 

3| ■ ■ ■ ■ — I 3| ■ ■ ■ ■ — I 3| ■ ■ ■ ■ — I 3| ■ ■ — 




-1 I ■ ■ ■ ■ — I -1 I ■ ■ ■ ■ — I -1 I ■ ■ ■ ■ — I -1 I ■ ■ ■ 

-40 -30 -20 -10 -40 -30 -20 -10 -40 -30 -20 -10 -40 -30 -20 -10 



Figure 18: Various shapes of EQ-inclusions when varying Im with b = 0.1 + 2i and Re 
6/a2 = -10. 
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